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Abstract 

This paper concerns the first passage times, denoted by r a , of Bessel processes to 
a point a > 0. We are interested in the case when the process starts at x > a and 
compute the densities of the distributions of r a to obtain the asymptotic forms of them 



as t — > oo that are valid uniformly for x > a. 



1 Introduction 

t^- ■ 

This paper concerns the first passage times, denoted by r , of Bessel processes of order 
v G R to a point a > 0. We are interested in the case when the process starts at x > a 
and compute the densities of the distributions of r a to obtain the exact asymptotic forms 
\ of them as t — > oo that are valid uniformly for x > a for all orders v. If v — ±1/2, we 
have well-known explicit expressions of them, which are often used in various circumstances, 
while otherwise there is quite restricted information on them that is available. In the case 
when < x < a the distribution of r a solves a boundary value problem of the associated 
second order differential equation on the finite interval (0, a) and the distribution of r a or its 
density is represented by means of the eigenfunction expansion (p], [I], [8], etc.). In the case 
x > a, however, the region for the differential equation is the infinite interval (0, oo) and the 
corresponding representation becomes a generalized Fourier integral, which it seems hard 
to derive an asymptotic form of the density directly from and there have been only a few 
partial results as given in [UJ, [13], where v — and relative ranges of x are restricted (in 
addition to the cases v = ±1/2). In a recent paper Hamana and Matumoto [5] have derived 
a certain explicit representation of the distribution of r a for the case x > a (as well as for 
the case < x < a) and computed an asymptotic form of the distribution and in [6] they 
have also given corresponding results for the density but the uniform estimate of the density 
does not readily follow from their result. For the present investigation we are motivated by 
the study of Wiener sausage of Brownian bridge in H d joining the origin to a point x £ H d 
over a time interval [0,t] when |x| grows linearly with t ([15]): the volume of the sausage 
swept by a ball of radius a is expressed by using the density for r a with x = |x| > a and 
v — (d — 2)/2. Although this concerns with the hitting time of mult i- dimensional Brownian 
motion only, it would be mathematically natural to study the problem in the framework of 
the Bessel processes. 



Let X t be the Bessel process of order v E R, whose infinitesimal generator is given 

by 

C <y) = + 2v + l . A. 

2 dx 2 2x dx 

If 2z/+2 is a positive integer, X t represents the radial part of the standard 2z/+2-dimensional 
Brownian motion. If v > — 1 we write <i for 2u + 2: 

d 

a = 2v + 2, or what is the same, v = — — 1; 

the process X t is sometimes called the <i-dimensional Bessel process. Let P x be the proba- 
bility law of the process X t started at x > and E x the expectation by P x . Let r a denote 
the first passage time of X t to a > and the density of the distribution of r a : 

q u (x,t;a) = j t P x [r a <t}. 

We also write for q u where d — 2u — 2, if v > — 1. 

In what follows we suppose > unless stated otherwise explicitly. In the end of this 
introduction we shall observe that there is a simple relation between q v and q~ u and the case 
v < is reduced to the case v > and vice versa. If u > 0, the origin is an entrance and 
non-exit boundary to the positive half line as is well-known. We shall use the two indices d 
and v interchangeably, understanding that they are related as above. Put 

p u t (x) =p ( t d \x) = (2-Kt)- d / 2 e- x2 ' 2t . 

For the process X started at the origin pf\x) is the density w.r.t. the invariant measure 
CdX d ~ 1 dx of the distribution of X t , where q = 2n d ^ 2 /T(^(d + 1)), the normalizing constant. 

In [13] the present author obtains the following result among others. Put k = 2e~ 27 , 
where 7 = — / °° e~ u lg u du (Euler's constant). 



Theorem 1 If u — 0, then uniformly for x > a, as t — > 00 
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for x 2 > t. 
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This theorem does not identify in any sense the asymptotic form of q^ 2 '(x,t; a) for x > 



y/ 4t lg lg t (see Lemma [TOl in Section 2). The objective of this paper is to complement this 
(Theorem HJ), and at the same time also to obtain an asymptotic form of q v for v > when 
the Bessel process is transient. For the latter we have the Green function, G^(x, y) say; we 
need to bring in G(x) = G^(x); = G^ d \x,0) = G^(0, x), or explicitly 



G{x) 



pf\x)dt 
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Theorem 2 Let v > 0. If v ^ 1, then uniformly for x > a, as t — > oo 
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In the case v — 1 £/ie same estimate holds true if the error term given by the O -symbol is 
replaced by 

P( 1 + I f /X2) ) fora< X <Vi; 0(!|i(^) 3/2 ) f„r X >V-t. 
Here (and here after) a V b = max{a, 6}, a A b = min{a, 6} for a, b real. 

Remark 1. (i) For random walks on the d-dimensional square lattice Z d we have analogues 
of Theorems [TJ and H] [12] . The form of the principal term in the formula is intrinsically 
the same as and in fact suggested by one for the walks. For d = 2 Theorem |2] (or Proposition 
[9]) provides an improvement of Theorem 1.4 of [12] in view of Theorem 1.5 of it. 

(ii) In Theorems [T] and [2] the phrase "uniformly for x > a, as t — > oo" may be replaced 
by more logically appropriate one "for x > a,t > 2". We have adopted the expression as 
given above to emphasize that the formulae asserted is meaningful only when t is large (the 
estimates are valid but poor as x — > oo with t being confined in a finite interval). The same 
remark may be applied to similar statements in which Landau's o symbol is not involved. 

(iii) In the proof of Theorem [2] we give a more precise expression of the error term, 
which shows that the order of magnitude cannot be improved at least for x < \rt. 



Define a function A u (y),y > by 

K(y) - 



(2tt 



v+l 



y>0 



and Aj,(0) = lii%j, A v (y). Here K u is the modified Bessel function of second kind of order 
v\ A„(0) is well defined, and in fact A(0) = and for v > 0, 



A„(0) 



2tt 



v+l 



G{1\ 



Theorem 3 Let v > 0. Then uniformly for x > as t — )■ oo, 



(d) 



(x,t;a)=a d - 2 A u (^y ( f\, 



x) l + o(l" 



(4) 



If v > 0, dl]) holds true as x A t — > oo while this is not true for v = 0, as is ascertained 
by referring to Theorems [1] and [2] 

The two-dimensional case is particularly interesting and deserves to state separately. 



Theorem 4 Uniformly for x > as t — >■ oo 

q( 2 \x, t, a) 



TV 



Pt \x) 



K (ax/t) 



(l + o(l)). 



Substitution from the formulae K (u) = — lg(ie 7 w) + 0{u 2 \gu) (u 1 0) and K (u) 



tx/Iu e~ u (l + 0{l/u))) makes the right-hand side above explicit if x/t goes to or oo. We 
shall actually compute errors in the formula of Theorem H] (Propositions |9] and ITS"]) so that 
for x > y/t (outside the parabolic region), 

q {2) (x,t,a) 7T / 1 \ 

= ^2%ax/te ax ^(l + 0(t/x)] if x/t oo. (5) 

The factor e ax ^ in the last formula may be understood to be natural by comparing with the 
Gaussian kernel pf\x — a). 

Theorem 5 For each v > 0, for t < 1 and x > a 

/ 1\ , , x — a ( (x - a) 2 \ /a\ ( d ~ 1 )/ 2 ?N 



qW(x ' i; a) = 71^ exp { - 2t) UJ (1 + 0{t/a »>• 

From the scaling property of Bessel processes it follows that 

q(x, t; a) = a~ 2 q(x/a,t/a 2 ; 1). 

For the proof of the foregoing theorems we shall mostly consider only the case a = 1 and 
write q(x,t) for q(x,t; 1). 

The estimation of q(x, t) is carried out in the following three cases 

(i) x < yt] (ii) yt < x < Mt (with M arbitrarily fixed) ; (iii) x/t — > oo, 

each being discussed in Sections 2, 3 and 4, respectively. The methods employed in these 
cases are different from one another. Roughly speaking, for the case (i) we apply the Laplace 
inversion formula and the Cauchy integral theorem, which provides somewhat finer estimates 
(Proposition [7]) than those given in Theorem [2j For the case (ii) we exploit the fact that 
any Bessel process of order v > — 1 can be decomposed as a sum of two independent Bessel 
processes and apply the result of the case (i). This gives some error estimate to the asserted 
asymptotic form of q in the case x/t — > (Propositions 191 and H2l . To include the case 
x/t — > v > an additional argument is employed. For (iii) we use the expression of q given 
in the following 

Lemma 6 Let (3 = (d + 1)(3 - d)/8 = {\ - u 2 )/2. Then 

(d-l)/2 



gW(M;a)=g) Eg" exp{/3/ o 



ds 



/o B 



T n = t 



yphxt? 



e -<*-ar/« ( 6 ) 



where the conditional expectation is taken w.r.t. the probability measure of the standard 
linear Brownian motion B t . 



Theorem [5] follows immediately from this lemma. If x/t is large enough, one can evaluate 
the conditional expectation in ()6]) for large t, of which the dependence on v comes only from 
P. Otherwise, however, a direct evaluation of it seems hard; our results rather give a precise 
estimate of it valid uniform in x. The proof of Lemma is given in the subsection 4.1. 

Throughout the paper C, C, C", etc. will be used to denote constants whose precise 
values are not important for the present purpose; the same letter may indicate different 
constants depending on the occasions where it occurs. 



We conclude this introduction by mentioning some simple facts for the case v < 0. The 
formula (j7]) given in the beginning of the next section entails that for v < 



.CM> 1 X 2M 



,2M 



(where the super script ^ designates the order of the Bessel process that is concerned) , so 
that 



q"(x, t; a) = g'"'(x, t; a) 



x 



2\u\ 



a 2 M ' 

The Bessel process of negative order visits the origin in a finite time with probability one 
and we have explicit formula 

q»(x,t;0) = A H (0)x 2 ^p[ l '\x) 

due to the integral representation K v (z) = 2~ 1 (z/2) u J* °° exp ( — u~ l — ^z 2 uju l/ ^ 1 du. By a 
comparison argument we have the inequality 

/ q v [x — a, s; 0)ds < / q u (x,s;a)ds (x > a, t > 0) if — 2 _1 < v < 
Jo Jo 

and the same one but in the opposite direction if v < —2~ 1 . 



2 Proof of Theorem [2 

For any v G R, 

^[exp{-Ar }] = ^r^f-u (A > 0, x > a > 0), (7) 

as is well-known and may be derived by solving the problem: C^'U = \U (x > a) with the 
lateral conditions U(a + 0) = 1 and U being positive and decreasing, of which solution is 
unique. (See (ITUj) below.) 

For v — 1/2 (i.e, d = 3) we have a particularly simple expression of q^: for x > a 



5 



ae -(z-a) 2 /2t , a 



^ 3 Hx,t;a) = ^ _ (l-~ 



which trivializes this special case of Theorems [2] and |3] and is helpful for making a guess at 
the asymptotic form of q in general cases. 

In what follows we let v > and, when there is no risk of confusion, we suppress the 
super-script ^ from q( d \x, t; a) and pf\x) except for the statement of theorems or lemmas. 

Put G x {x) = J °° Pt {x)e- Xt dt. We know 

G A (x) = ^ F ,(^j K u (xV2X) (9) 
([3], p. 146). It is convenient (and natural) to write the representation ([7]) in the form 

£„[exp{-AT 8 }] = (10) 



For v > let G(x) = lim^o G\(x), so that 



G(x) 



p t (x)dt 



r(f-i) 

2rr d/2 x d-2 ■ 



If v is not integer, K v {z) = [tt/2 sin(7rz/)J l~ v {z) — I v {z) , where l v is the modified Bessel 
function of the first kind of order v and given by 



£ 

k=0 



(z/2) 



2k 



T(u + k + l)T(k + l) 



'ir. 



for | arg^l < tc ([9J, p. 108). 

Proposition 7 Lei v > and M > 0. For a < x < M\/t, 



q {d) (x,t; a) 



G(a) 



d-2 



rf°(a) 



+ a 2 r]( y x/a, t/a 2 ), 



(12) 



wift t/ie function r](x,t), x > l,t > 2 admitting the estimate 



T](x, t) = O 



1 — X 

t»+ 2 



if v > 1; 



1 - X" 1 



t 2u+l 



if < i/ < 1 



and n(x, i) = O 



x 



-i 



t «/+2 



1st 



i/ i/ = 1. 



fSee / TTgj) . p7P and |H|) /or more exact forms ofrj(x,t).) 

Remark 2. One might suspect that the function q*(x,t) := [G(a)] _1 p|^(x — a)[l — 
(a/x) d ~ 2 ), an analogue to the exact form of q^ 3 \ can take place of the leading term in the 
formula ( I12p . Since the difference of them is at most the magnitude of 0(x^ 2 ~ 2u ^ yl /t u+2 ), 
this is true if v < 1; in the case v > 1, however, the difference becomes much larger than 
r](t, x) as x gets large, so that the replacement causes a larger error term. 

Proof of Proposition [7| The Laplace inversion of ffTUj) gives 

3 GJx) 



q(x, t; a) 



1 

27ri 



GJa) 



e tz dz. 



(13) 



For the evaluation of this integral we follow the argument made in [12] for the random walk 
of dimensions d > 3. Motivated by it we decompose 



G x {x) _ G x {x) - (a/x) 2u G x (a) a 2u 



Gx(a) 



G\(a) ./• 
G A (x) - {a/x) 2 »G x {a) 



2v 



G(a) 



R(X;x) 



(14) 



where 



R{X;x) 



1 



1 



G x (a) G(a) 



G x (x) - (a/x) 2iy G A (a) 



;2v 



+ 



2;/ 



By (|9]) the contribution to (|T3|) of the first term on the right-hand side of (IT4"|) equals 

1 



G(a) 



Pt(ar) - (a/x) v pt{a) 



The error term in Proposition [7] is then written as 



[ rioo 

rj(x,t) = — R(z,x)e tz dz. 
2m J-ioo 



(15) 



As z — > oo 

K v {z) = (7r/2z) 1 / 2 e -*(l + 0(1/ z)) (| arg^| < tt). (16) 

Hence G z (x)/G z (a) = 0(e~( x ~ a ^^) as z — > oo (with x fixed) and we see that the function 
R(z, x) rapidly approaches zero as z — > oo in the sector | argz| > tt — 5 for any 5 > 0. This 
together with the fact that K u (z) has no zeros on the right half plane 3ftz > (cf. [16] : p. 
511) (hence G z (x) has no zeros on — tt < aigz < tt) permits us to apply Cauchy's integral 
formula to transform (ITS"]) into 



Vfat) = — 

2m Jo 



R(—u + iO, x) + R(—u — iO, x) 



e~ tu du. 



(17) 



Let v be not an integer. We show that uniformly for a < x < M\fi, 



r](x,t) 



-v-\ 



2»+ l (u- 1)I» 



1 x 
1 - — + - 



Br 



1 — X 



-2v 



with 



\r(t, x)\ < 



C 



t 2u+l 



X 2 \ 1 

— V — 

t J t v 



X 



- t v+2 t \ X 



T(-2v)t 2 »+ l 

(0 < v < 1) 
(i/>l), 



- I + r(t,x) (18) 



where the constant Bq is given in ( l20j) below and l/r(— n) = if n is a non-negative 
integer. The estimation 77 is simple apart from the uniformity in x, which we must take 
care of in dealing with the dependence on x. Let a = 1 for simplicity. Recall the definitions 
of G\(l) and G(l). From the power series expansion of I„ given in ( ITT]) we then deduce 
straightforwardly 



G(l) 



2A)^[/_ i ,(v / 2A)-4(v / 2A)] 



lim A40 (v / 2A)-/_,(v / 2A) 

r(1 -"(£ r ( J^i) M 



+ 



2^ ~ (A/2) fc 

&r(i/ + * + i)fl 



l + AiA + A 2 A 2 + 



A !/ 



B + SiA + ---] 



with 



Ax = [2(1 - v)}-\ B = 2-T(l - i/)/r(z/ + 1). 
Let /(A) = G A (1)/G(1) - 1. Then 



G A (1) G(l) G(l) 



1 + /(A) 



/(A) + [/(A)] 2 



G(l) G A (1) 



(19) 
(20) 

(21) 



Also, noting 



a^Gx(x) = G 



we obtain 



x 2 »x 



G x (x) - (l/x) 2 ^G A (l) G, 2A (1) - G A (1) 



G(l) 



G(l) 

A x A(x 2 - 1) - B Q X u (x 2v - 1) + H\(x), (22) 



where H\(x) 



A 2 X 2 (x 4 - 1) + • • ■ 1 + [fixA^+^x 2 ^ 2 - 1) + 



the remainder term. Since 



/(z) = Od^l 1 '/ 2 x / 4 ) as 2 — > oo along the negative real line, the equality ( 1T9|) entails 

/(A) = AiA - 5 A" + 0(|A| (l//U)+1 ) (|A|<1) and ± i0)| < C\u\ u ' 2 (w > 0). 

From the power series expansions of I± v (z) it follows that A k = 0(l/T(—i / + k + l)k\), B k 
0{l/T{u + k + l)k\). 

With these preliminary discussions we now compute the integral (IT7|) . First consider 

11(A) := [/(A)] 2 [Ga;2 A (l) — G A (1)]/G A (1) 
and observe that for u > 0, 



+ iO) - Ti(-« - i0)| < G u 2{l/A1) u u (x v — 1) + Gin^ + (i 



,2iAf2 



!) + •• 



+ G w 2l/ «(x 2 - 1) + C[u 2 (x* - 1) + 



with certain constants C k ,C' k that are dominated by a constant multiple of 2 k k 3u /(k\) 2 . 
Here we have exploited the fact that the terms of integral powers c n X n involved in 7i(A( 
cancels out in the difference on the left-hand side. Employing the simple inequality x s — 1 < 
(1 V s)(l — x~ 1 )x s valid for all x > 1, s > we infer that for 1 < x < M\/t, 



/ TU-u + iO) -Ti(-iz + i0) 
Jo 



< C 1 



<C"(l-i 

x 



£2(VA1)+1 



fc=0 



t u+k k\ 



Y °° k 5u+1 4 k x 2k ^ 
5^ 



fe=l 



x' 

— ) + 



t 2 ^ 1 V t 



x 



£2(wu)+i y t 

where G' depends on M. Next, in the same way we see that if T 2 (A) = f(X)H(X), 

l\/x 2 \ u + 1 1 



(23) 



|T 2 (-u + z0) - T 2 (-u - iO)\e- tu du < c(l - ^ 



^Al+l 



(24) 



Thirdly let 



F(X) = -f(X)(A 1 X(x 2 - 1) - B A V - 1)), 
which in view of f ll4p . §2U§ and (T52J) would involve the principal part of -R(A; x)x 2v . Then, 

F(X) = AiBoX^ 1 ^ 2 " - 1 + x 2 - l) (1 + dX + ■ ■ •) - BlX 2v [x 2v - + C[X + ■ ■ •) 

apart from the difference of an entire function, and for s > — 1, 

1 



1 

27ri Jo 



- (-u + i0) s + (-u-iOy 



e~ tu du 



V{-s)t 



s+1 ' 



(Here the identity T(l + s) siri7rs = — 7r/T(— s) is used; remember that l/T(—n) = if n is 
non- negative integer.) Hence 



— / \-F(-u + iO) + F(-u - iO)}e- tu du 
2m Jo 



A l B { 
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2u 



X 



1 



Y(-v- l)t"+ 2 



x 2v - 1 



l + 0(l/t) 



(25) 



Finally collecting the bounds (TIT]) and f fT9|) through ff24|) (of which we divide each formula 
by x 2u since we have multiplied it in f )22|) ). noting B /T(— v — 1) = {y + l)/2 v T(v) and 
making elementary comparison of terms that appear on the right-hand sides of them we 
find (HBl) to be true. 



Let v be a positive integer. The arguments are similar to the above. In place of ( fl9l) we 
have 



Ga(1) 
G(l) 
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~o~ ) - hf 7 + Wo E 



(A/2)' 



fc=0 



fe!(i/ + ife)! 



+ 



0(A) 



= "l + AxA + A 2 A 2 + ---] - A i/ (lgA)[5 + 5 1 A + -- 
where g(A) = ao + aiA + ■ ■ • is a certain entire function (with do = 1 for v — 1), and 

(-ir 



" 2 " [1 + 1 ? ( r 27)] and 5o 



-[2(1/ - 1)]- 



(i/>2) 



2"z/!(f - 1)!' 



(26) 



By noting 



(2tt)- 1 3 - [lg(-M + 20)] fe + [lg(-u - «0)f = -fc(lgw 



(fc = 1,2,« > 0) 



we compute the integral in ifPTj) to see that if ^ is an integer greater than 1, then 



r](x,t) 



1)^^ + 1)1^0(1- — 



x 



1 \ 1 



0\ 



'1 - x- 1 
t»+ 2 



X 



(27) 



(The coefficient of the leading term coincides with one appearing in (fTBl) . so that ( TT8|) 

. . ..) \i v — \ and / is defined as before, the leading term of 
/G(l) being (5 2 A 2 lg A)(-x 2 lg(x 2 A) + lg A), we see 



is valid for f 



-/(A) 



2,3,4, 



r?(x,t) = -45 — + O) 



1 — x 
t»+ 2 



(2f 



The foregoing two formulae obviously implies the desired bound for r](x,t). The proof of 
Proposition [7] is complete. □ 



Proposition 8 Let v > 0. T/ien the function r)(t,x) defined via ITU^) admits the estimates 

and r/(t,x) =of^| (-^) ^ if v = 1. 
i/iai are i>a/z<i uniformly for x > \/t > 2 . 



Proof. We can proceed as in the preceding proof except that in place of (1221 we make 
decomposition 



■2v 



X 



(G x (x) - (l/x) 2 »G x (l)) = G X 2 X (1) - G(l) + [G(l) - G A (1)] 



and estimate the contributions of the three terms on the right hand side separately. Let v ^ 
1. It follows from the preceding proof that the contributions of the last two terms to rj(x, t) 
is bounded by a constant multiple of x~ 2v t~ v ~ x = t~ 2v ~ x {\fi jx) 2v . As for the first term, 



on the one hand, we recall that G 2 (l) = C v {2z) v l 2 K v {^2z) = G^-^e-^l + o(l)) as 
z — > oo to see that 

^- 2 -|/(-^ ± iO)G-_^ ±i0 (l)| < Cu Xh »x- 2v {x 2 uf v -^^ = Cu 1 * v+V (a?u)-''' 2 - 1 ' 4 

for u > 1/x 2 (/ is the same as before) and the integration over u > 1/x 2 of e~ tu times the last 
member yields the magnitude of the order 0{t~ <ylhu ^ +u+l (\ft/ x) v+l l 2 ^j . On the other hand 
we have x~ 2u \f(— u ± iO)G_ x 2 u±i0 (l)\ < Cu 1Au x~ 2u for < u < 1/x 2 and the corresponding 
integral does not exceed the foregoing magnitude. By ( ITT)) we find the asserted bound for 
v ^ 1. The case v — 1 is omitted, it being similarly dealt with 



3 Proof of Theorem [3] when x jt is bounded 

In this section we prove the formula of Theorem in the case when xjt is bounded. It is 
convenient to treat the cases v — and v > separately. In both cases one may suppose 
xjt — > v and the subcases v = and t> > are also separately treated since different 
arguments are employed for them, although the framework is the same. In the case v = 
we shall provide estimates of the error terms, that are not given in Theorem [3j 

3.1. The case v = 

In the course of proof of Theorem [3] we shall prove the following proposition, which 
entails the formula of Theorem [3] in the case xjt — > 0. 

Proposition 9 Let v = 0. It holds that uniformly for \ft < x < t/2, as t — > oo 

(2) / \ \ 

m 0*0 \ 



<1 



(2) (x,t,a) 
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K (ax/t) 



pf\x)+0\ 



\\g{t/x)f + l t 



For the proof of Proposition we shall use Theorem [JJ which it is convenient to reduce 
to the following slightly weaker form. 



Lemma 10 Let v = 0. Uniformly for x > a, as t — ^ oo 



q^ 2 \x, t; a) 



2ttp?\x) 



2np i2 \x) 



lg((z/a) 2 ) 
{\g{t/o?)Y 

lg((x/a) 2 ) 



+ 



(lg(t/a 2 )) 2+ \lgt 



Ost) 
i 



for x 2 < 2tlg(lg£), 
for 2t lg(lg t) < x 2 < At lg(lg t) . 



In what follows we let a = 1. We need the following lemma from [T3] in which v = 0. 
Lemma 11 For any v, there is a constant c = c u > such that for all x > 1 and t > 1, 

q {d \x,t)<cp^ 1 (x). (29) 



Proof. The proof follows from the parabolic Harnack inequality as in the case v = 0. 

We use the fact that the Bessel process of order v = is the radial process of the standard 
two dimensional Brownian motion B^ . Let x denote a generic point of R 2 and M ( 2 ) the 
probability of B^ started at x . We can suppose that the initial point B^ = x is on 
the upper vertical axis so that x = (0,x). Write £ t and Y t for the horizontal and vertical 
components of B\ , respectively, and let Tk be the first hitting time of the vertical level K 

(2) 

by B\ ': T K = inf{t > : Y t = K}. Then the space-time distribution of (T K , £ Tk ) is given 
by 

P^T[r K edt^ K eam x-K m (1) 

dtdt = — Pt {x ~ K)pt ® (30) 

(cf. [7J, page 25), which yields the representation 

q(x,t) = fds r J? + (x~ K) 2 )qUe + K\ s)d£. (31) 

K may be any number between 1 and x = |x| but we suppose 4 < K < x/2. With a fixed 
K we are to compute the repeated integral on the right-hand side by using the formula of 
Lemma [TUl It is remarked that we make no use of Lemma [TU1 in the case when x/t — Y v > 0. 

We break the rest of the proof into four parts. For the case x/t — > certain elaborate 
computations directly yield the desired formula of Proposition [9] (Parts 1 and 2). In the case 
x/t —7- v 7^ 0, on the other hand, we first show the existence of limit of the ratio q(x, t) /pt(x) 
(Part 3). While it is difficult to identify the limit along the same line as in the case x/t — > 0, 
with its existence at hand another way determines the limit (Part 4). 

Throughout the proof we suppose that for some M > 0, 

y/i < x < Mt. 

The constant K > 4 may be fixed arbitrarily prior to Part 3, in which we need to take K 
large enough, so we do not assign K a specific value. We write p t , q for pf\ q^ to be 
consistent to our convention that the super-script ^ is dropped, while we continue to write 
Pt^ if k 7^ d. We put 

J b J — oo l S 

the contribution to the integral in (|3T|) from the interval b < s < c. 

Part 1: Estimation ofI Ctt - Here c is a constant not less than 4. In the identity pa 2 +qj3 2 = 
pq(a — (3) 2 + (pa + q(3) 2 , where p,g 6 R with p + q = 1 and a, (3 may be vectors in any 
Euclidian space, take p = (t — s)/t and divide both sides of it by 

T = pqt= S -^-^ 



to obtain 

\o? + ^/ = \{P-a) 2 + ^(a+ S 1 {P-a)) . (32) 



Then substituting the two-dimensional vectors a — (£, K), ft — (£, K — x) leads to 
Pt-s(^ 2 + (x-Ky)p s (^e + K 2 ) = ^Pt{x)e-?l 2T e* V \-^(K-- t x^ 



(33) 



In the repeated integral of I Ci t we split the range of integration w.r.t. £ at £ = ±\/4slg lgs. 
We claim that 

let = Pt{x) / ~, t—!= exp < — —= \K — -x 



x 



(t-s)y/T \ 2TV t 

/■ ig(£ 2 + ^ 2 ) e _? /zr 

Jy/e+K 2 <y/4sl S l S s (\g S ) 2 Vr 

1 



+ i2(s,*)+0 



(lgs) 



ds, (34) 



where _R(s, t), the term that comes from the remainder term in Lemma fTUl is o(l/ lg s). For 
the part |£| > ^/4slgs we have only to substitute the expression of q given in Lemma [TOl 
and apply fl33l) (note that the bound in Lemma [TOl actually holds uniformly for t > 4 simply 
because q(x,t) is bounded there). For the integral on |£| > V^iTgs we need to take p s+ i in 
place of p s in ( 1331) so that the corresponding contribution to J Cjt becomes 

Pl+lW i' (^7F exp { " 2^(* " ITT 1 ) 2 ) x 

e -5 2 /2T' 

-d^ds, (35) 



where T' = (s + 1)(£ — s)/(t+ 1); since the inner integral is 0(l/(lgs) 2 ) uniformly in t, a 
simple change of variable gives the error term in (|34[) . Thus we have verified the claim. 

Scaling the variable £ by y/s and dominating 'l 2T by t~^l 2s if necessary we see that 
the quantity in the big square brackets is evaluated to be 

v27T o( } )+R(s,t). (36) 



lgs V(ig S )2^ (t _ s )/ t 

We must compute the integral 

x-K f 1 / s \ 2 ) v 7 ^ , 

J : = / =■ exp < 7z\K x) K ds. 37 

Jc (t-s)VT P 1 2T\ t J J lgs 1 ; 

Now we suppose x/t < 1/2 so that Igt/x > lg 2 and take c = 4. By a simple change of the 
variables of integration, e.g., according to u — (x/t)y/s, which transforms s/t to u 2 t/x 2 one 
can easily find that J ~ p t (x)ir/ \g(t/x) as a?/t — >■ 0, t/x 2 — > (which is enough for Theorem 
[3] restricted to the case x/t — > 0). But this way does not give the error estimate asserted 
in Proposition |9j To improve the evaluation of the integral we transform the variable of 
integration by 

s ( to \ , . 

o = - i.e., 8 = —-\ 38 

t — s V 1 + a / 

entailing the relation ds = (t — s) 2 dcr/t = (t — s)^Jt /ta da. 
Noting the inequalities 

11 2 1 sx x 2x (2 x 2 a 

< — < and = < — V 

T ta t — s T t t — s t \x t 

we may write the exponent appearing in the integral of ( 13?]) in the form 



where 5 = 5(t, x, s, K) is a function of t, x, s, K that satisfies < 5 < AK/x (provided that 



K < x). Further transform the variable a to u = xya/t. Then 

xds xda 



(t 



T 



to 



2du, 



and we obtain 
J 



2V27T 



(x/t)y/c/(i-c/t) lg (t/x) 2 m(u) 



■ exp 



2^2 



K 2 x 
2t 2 u 2 



u 



1 -5)— )du 1 + - 



.r 



(39) 



where m(u) = u 2 / (1 + u 2 t/x 2 ). We apply the inequality |1 — 1/(1 + r)| < \r\ + r 2 / (1 + r) 
(r > —1) with r = [lgm(w)]/ lg(t 2 /x 2 ) for which (1 + r) _1 < (lgc) -1 \g(t 2 /x 2 ) to see that 

2V2^ V2^ / ! + [\gm(u)} 2 \ 

\g[(t/x) 2 m(u)}~ \g(t/x) + \ [\g(t/x)} 2 ) 

uniformly valid if u is confined to the range of integration. Owing to the identity 

e~ b / 2u2 - Xu2 l 2 du = t/^e"^ (A > 0, b > 0) (40) 

Jo V 2a 

and the bound / °° | \gm(u)\ 2 e~ u2 ^ 2 du < C the formula (1391) reduces to 

J=[7r/lg(t/x)][l + 0(l/lg(t/a;))]. (41) 

Taking the computation carried out right above into account one also observes that the 
contribution of the error term in (136|) is O (p t (x) / [\g(t / x)] 2 ^ and concludes that uniformly 
for y/t < x < t/2, as t — )■ oo 



let =Pt{x) 
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lg(t/x) 



1 + 



\g(t/x) 



R 



(42) 



where R, the term corresponding to R(s,t), is o(l). 

Part 2: Estimation of I Q c . The integral Io jC (t,x) is dominated by 

r < c ] SU P ^ Z ^ft-.(>/€ 3 +(a:-iO a )de 

J-oo V? 0<s<c t— S vv y 



< Cy/-ce- K ' 2 l 2c ^ Pt {x-K) < [Ce-^^e^ 



Pt(x)-, 



where for the first inequality we have used the bound P y \j\ < c) < C\fce~( y ~ 1 ^' 1 l 2c / (y — 1) for 
y > 1, a bound obtained from the one-dimensional result (cf. Lemma 3.2 of [H]). Combined 
with (T42D this shows that q(x,t)/p t {x) = 0(1/ Igt) at least for y/2t lg lg t < x < y/4t lg lg t. 
Using this bound instead of the second one of Lemma [10] we obtain 

R(s,t) = 0(l/(\gsf) 

so that R = 0(1/ \g(t/x)). The proof of Proposition [9] is complete. 



Part 3: Proof of convergence. Here we suppose x/t — > v > and prove that there exists 
lim q(x,t)/p t (x), the limit value depends only on v and the convergence is locally uniform 
in v . Here we use Lemma [TT1 but does not Lemma [TUl With the help of fl33|) it gives 



where s' = s+1, t! = t+1 and (1331) is applied with p s+ i (V^ 2 + K 2 J in place of p s (\/£ 2 + .fT 2 
(see fl35]) ). One observes that the integral above is at most 0(e~ vK ^) (use e.g. (HQ])). Also 
a quite crude estimation shows I t / 2>t (x,t) < Cpt{x)e~^ 2K ~ xS> / 8 *. Combined with the result 
of Part 2 these show that for any e > one can choose K large so that 



lim sup 

t— >oo,x/t— >v 



q(x,t) - I c>K 2(x,t) 



Pt(x) 



< 6. 



(44) 



Define s) by 

hK(Z,a) = q(y/? + K*,s)/p a (y/e 

By Lemma ITT1 

Iik(£, s) < C exp 
and keeping this bound in mind we see that 

Ic,K*{x,t) 



K 2 



2 S (s + r 



Pt[x) 
1 

p t (x) 



K 2 



ds 



°° x - K 



t 



p t ^ s + (x - k) 2 ) Ps ( + k 2 ) h K (t, S )dz 



R2 ds r°° x- K 



2nT 



t 



exp 



K 



X \2 



k 2 ds 



2ns J - 



cxp 



2TV- t 

(k-vs) 2 e 

2s 



hK{£,s)d£ 



2s 



hK{£,s)d£ 



as x/t — > v. This together with (j4*4"l) shows that q(x,t)/p t (x) is convergent and the limit 
value does not depend on the manner of x/t approaching to v. The required uniformity of 
the convergence is easily ascertained from the arguments made above. 

Part Jf.: Identification of the limit. The proof rests on the identity 



p t {x) 



q(x, t — s)p s (l)ds. 



(45) 



which follows from the Markov property of the Bessel process and also from the identity 
ffTUj) . We may suppose that x — tv, v ^ 0. By Part 3 

q(tv,t - s) = Xpt- S (tv)(l + o(l)) as s/t -> Q,t -> oo (46) 

for some constant A = X(v) > 0. Since 



Pts{tv)p s (l) 



1 



Substitution of 



p t (tv) 2vrs(l - s/t) 

into (141)1) yields 

1 



exp 



u 2 s 



2(1 - s/t) 2s 



1 

A 



lim 



1 



p t (tv) Jo 



2vr Jo 6XP 



p t _ a (to))p s (l)ds 

v 2 s 1 \ds 
' ~2~ ~ ~2s)~ 



(47) 



(see [3], Eq(29) in page 146 for the last equality). Hence A = ir/K (v) as desired. This 
completes the proof of Theorem [3] in the case when x/t is bounded. □ 



Remark 3. In the case x/t — > v > 0, to compute the value lim q(x,t)/p t (x) along the 
same line as in Part 1 is hard since our knowledge of the behavior of q((C,,K),s) is poor 
for small values of s that significantly contributes to the integral of (l3Tj) . This point would 
be well understood from the argument made in Part 3 above. One notes that from Part 
2 we know only that Io^ c (x,t) becomes small if K/c is large enough, while the error term 
0(l/(lgT) 2 ) in the estimate (134]) depends on c . 

3.2. The case v > 0. 

Proposition 12 Let v > 0. It holds that uniformly for \/t < x < t/2, as t — > oo 
qW(x,t,a) 1 r^(( x \ d ~ 2 



P <-»W G(a) - V«^>J <°<" <!/*). 

= 0(x/t) (u > 1/2). 

We use the fact that the square of Bessel process of dimension d > 2 is the sum of those 
of two independent Bessel processes of dimensions 1 and d! — d — 1 ([ID])- Let (Yt) be the 
one-dimensional standard Brownian motion started at x and (£ t ) the Bessel process of the 
dimension d' started at and independent of (Y t ) . Then the law of the process (X 2 ) is the 
same as the law of (£f + Y 2 ). Let T% = inf{t > : Y t = K}. Then in place of ( l30l) we have 

P\T K e dUr K e d(] = *_K p?[x _ K)p ^ i0cd ,^ 

dtdt, t 

so that 

q(x,t) = fds r ^— -p?2 s (x - K)ptl{i)qUe + K\ s)c d ^ d - 2 dC (48) 

JO J-oo t — S K v ' 

The proof of Proposition [12] given below is analogous to the one given for v = and we 
proceed parallel to the lines of the preceding proof. 

Part 1: Estimation of l c ±. We write Ib, c (x,t) as before for the integral in (T48l) restricted 
on the interval [b,c]. The product Pp2 s (x ~ K)pi-l{0 appearing in the integrand may be 
written as 



Pt-s^e + {x-KY)= p t ^ s {x - K)e-?W- S \ 
which we further rewrite in the form 



1 ^ d/2 p t (x - K) exp ( - {X ~ K ? s ) e -ew-). (49) 



t-sj I 2t(t - s 



We split the range of ^-integration at £ = ±a/4s lg s in the repeated integral of I C)t - 
The integral on |£| > ^/4s\gs is disposed of by employing Lemma fTTl as before (see (1351) ). 
As for the integral on the other part we first evaluate the contribution of the term (£ 2 + 
K 2 )~ u p s {l) /G(l), which, on using f T49|) . is dominated by a constant multiple of 

TO ' L («-s) J / 2+1 F l 2t(t-s) I 4i<VITii7 



It is convenient to split the outer integral at s = t/2 and let R± and Rf be the parts 
corresponding to s > t/2 and s < t/2, respectively. By performing the ^-integration and 
changing the variable by u = t — s we deduce 



/■t/2 

< Cp t (x) 

Jo 



X 



e~ x2/Au du < C"p t (x)x 2 - d . 



For the evaluation of Rf we replace the integrand by unity in the inner integral and have 
the bound 



Rf < 



Cp t (x)x r 1 / 2 A/lg s 



t 



exp < — 



2 

X s 



ds. 



s (d-l)/2 

Since the integral on the right-hand side is evaluated to be 0(x/t) d ~ 3 ^Jlgt/x^J or 0(1) 
according as v < 1/2 or v > 1/2, by taking account of the estimate for R> obtained above 
we deduce that 



x\ 



d-2 



R 1 <C'p t {x)(-) Jig- ifi/<§ and R x < C'p t (x)- if v > \. (50) 



tj 



x 



Let < v < 1/2. Then, in a similar way to the above, we evaluate the contribution of 
the error term in ([3]), denoted by R2 and make decomposition R2 = R£ + R£ ■ F° r R2 we 
note that / R e- ?2 / 2 (*- s ^ d - 2 ^ = 0((t - s)^ 1 )/ 2 ) and deduce that 

r t/2 

iRol < Cp t {x)xr v / u- 3/2 e~ x /Au du < C'p t {x)t- p - 
Jo 

also 



1^1 < 



Cp t (x)x r*/ 2 e- {x2 ' 2t2 >ds rV^ s 



J c s d/2 + u J Z d $ 



Cp t (x)x r l (lgs) 



t (U.^l)/2 



c S 



(d-l)/2 



e -^/2t*)s ds < C - 



d-2 



so that 



|# 2 i<C> t (z) - Wig-. 



(51) 



(52) 



Putting T = s{t — s)/t we have in place of (1331) 
^(A/e + ^-^kfv/e + A' 2 ) =p t (x)p¥\0 



V2nT 



exp 



2T V t 



Applying this together with ( 150]) and (1521) and making use of Lemma [TT] in the same manner 
as before we find that 



p t {x) 



l c,t 



x-K 



G{1) Jc (t-s)V2^T 



exp 



2T\ t 



x 



#(e)^e d - 2 ^+o(- 



d-2. 



y- v /s 5 +^ i< v^ 



+ 0^)7 \ fe- • 



d-2 



The quantity in the big square brackets may be evaluated to be 1 + 0(l/s[(t — s)/t]^ _1 ^ 2 
In order to evaluate the whole integral we employ the transformation fl38l) and follow the 
succeeding arguments up to flMj). We then conclude that 



Pt{x) 
G(l) 



x 



Let v > 1/2. Then, the integral of the third member in f l5T|) becomes bounded, so that 
we have [R^l < C'p t (x)x/t in place of the bound given therein. The other computations 
may be carried out in a similar way and we obtain I ct = \p t (x)/G(l)](l + 0(x/t)). 

Part 2: Estimation of Jo,c- The same computation as before gives the same bound of 
Jo,c (but here p t 
has been proved. 



Pt ), which is sufficient for the present purpose. Thus Proposition 



Part 3: Proof of convergence. The proof is quite similar to the one given for v = 0. The 
bound (143]) and hence the relation (jHJ) holds true without any alteration except that here 
q and p t are defined with d > 2. Define hx{C,, s) as before. Then 



Ic,K*{x,t) 
p t {x) 



K 



ds 



{2ns) d / 2 



exp 



(k - vsf e 



2s 



2s 



h K (^s)c d ,Z d - 2 dt 



as x/t — > v and as before we conclude the desired convergence. 

Part 4- Let Xd{v) be the limit of q(x,t)/pt(x) as x/t — > v > 0. The functional equation 
(J32J) holds true for all v > in view of the corresponding identity for the Laplace transforms. 
In place of ( H7|) we then have that if x/t — > v, then 



lim 



Pt- S {tv)p s (l) 



ds 



X d (v) Jo Ptitv) 

so that Xdiy) = A d (v). Thus we conclude the formula of Theorem [31 



v 2 s 



ds 



(2n) d / 2 Jo ^\ 2 2sJs d / 2 
2v d/2 - 1 K d/2 _ 1 (v)/(2nf 2 , 



□ 



4 Proof of Theorem [3] (the case x/t — > oo) 



In the following result we slightly refine the estimate asserted in Theorem [3] in the case 

x/t —7- OO. 



Proposition 13 For each v > 0, as both x/t and t tend to infinity 



qW(x,t) 



x 



V2rt* 



exp 



2t 



l + 0[- 

x 



(53) 



Our proof rests on the formula ([6]) and is broken into three parts. We begin with verification 
of the formula (Q. In view of the formula ([6]) the situation is distinguished whether v is 
less than 1/2 or not. The cases v > 1/2 and v < 1/2 are treated in the second and third 
parts, respectively. 

4.1. Proof of Lemma H 



Let 



j(x) 



d- 1 

2x 



By the Cameron-Martin-Girsanov formula 



f g (d) (x, s)ds = P x [t-h<r a <t} = E BM [Z(r a ); t - h < r a < t] 

Jt-h 



(54) 



for < h < t, where E BM indicates the expectation w.r.t. the probability measure of the 
one- dimensional Brownian motion B t started at x, r a is the first passage time for B t and 



By Ito's formula 



Hence 



and, noting 



L 



Z(t) = eIo^ Bs)dB ^^o^ {Bs)ds . 



' 1 m B t 1 ft ds 



Z(r n 



a \( rf - 1 )/ 2 



cxp 



(t < T ). 

(d-l)(3-d) [ T « ds 



o Bl 



P BM [t -h<r a <t} = f q {1 \x - a, s)ds 

Jt-h 

we deduce from fl54l the equality 

■ a \{d-l)/2 



qW(x,t) 



x 



E 



T~a = t q^\x — a, t) 



(p = (d- 1)(3 - d)/8 = (1 - 4z/ 2 )/8), which is the same as 

4.2. The case v > 1/2. 

Let v > 1/2 so that < 0. It suffices to show 



□ 



BM 



B f ' B~ 2 ds 

e Jo s 



T a = t 



1 + 0(t/x) 



(55) 



as x/t — > oo. Put 

v = x/t and J t (z) = £f M 
Then, by the strong Markov property of Brownian motion 



B f* B~ 2 ds 

e Jo s 



T a = t 



Ux) = / J 1/v (y)E BM 

J a 

By bringing in the measure 



jj;- 1/v B-*d S . 



i 



;r a > t ,B t _u v E dy 

v 



q {1) (y - a, i/v) 

— a, t) 



/i(dy) = HtA d y) 



q w (y - a, i/v) 

g( 1 )(x — a, t) 



P x BM [Bt-i /v edy], 



this may be written as 



J Ax) 



Ji/v(y)E BM 



rt-l/v 2 1 

e^Jo ;r a >t 



t-i/jj 



a; 



(56) 



An elementary (but careful) computation shows that 



fi(dy) 



y-g 



exp 



\[y-(a+ 1)] \ - |(1 + o(l)) ) (1 + o(l))dy (57) 



with o(l) — > as v — > oo uniformly in y > a, entailing that fi converges to the unit measure 
concentrated at y — a + 1 in the limit as v — > oo. 

For each non-random t > the conditional law P^ M [- \B to = x] is the same as the law 
under P G BM of 

(6(1 - s/to) + s(x - B t0 )/t + B s ) Q<s<to . 
Combining this with the well known fact that may read 

P™[B t + vs<a + 2~ 1 vs for some s > 0] = e~ v , 

we infer that 



(5f 



E, 



BM 
a+1 



1 1 1 

B s + v s > a + 2~ 1 vs for < s < t ;r a >t B t _ 1/v = x =l + 0(l/v). 

V V 1 J 



If the event in this conditional probability occurs, we have J * B s 2 ds < 2/av. Hence 



E. 



BM 
a+1 



< ^ Bs ds ;r a >t-- 

v 



Bt-i/i 



x 



l + 0(l/v). 



This remains true if the initial position a + 1 is replaced by any y > 1 + a/2. Therefore, 
with the help of (1571) and the trivial relation J\/ V {y) = 1 + 0(l/v) valid uniformly in y > a 
we conclude fl55l) from ([56]) • 



4.3. The case < v < 1/2. 

Here we have /3 > and we must evaluate the conditional expectation appearing in 
(1551) from above, the lower bound being trivial. To this end we apply the Kac formula and 
resort to the exact solution of a certain differential equation. In view of f|56|) and the result 
mentioned right after fl57j) it suffices to show that for all y > a 



W := E. 



BM 



exp <^ (3 



ds 



T a > t 



B 



t-l/v 



X 



< 1 



Co 

V 



provided v is large enough, where v = x/t as in 4.2. 

In order to obtain a tractable upper bound of W we discard the condition for non- 
absorption up to time t — 1/v and at the same time replace B s by (B s V a) in the integral 
in the exponent: also, we express the conditional expectation by means of the uncondi- 
tional realization of Brownian bridge given in ( 158]) and thereafter restrict the range of the 
expectation to the event 

\B t -i/ v \ < \fv t. 
Then, using the monotonicity of the function x V a we obtain 

ds 



W < E* M 



exp W 



(B s + (v- 2^)s) V a 
for all sufficiently large v. For a positive number t>* put 

ds 



+ Ce W/a 2 )t-vt/2 



U(y,t;v* 



E 



BM 



eXP Wo' 



[(Bs + v*s) V a}^ 



(t>0,yER). 



Then W < U(y, t; v*) + Ce~ x ^ if v* < v — 2y/v and we have only to show that U(y, t; v*) = 
1 + 0(1 /v*) uniformly for y > a. 

For each fixed, the function U(y,t) = U(y,t;v*) is a unique solution of the parabolic 
equation 

i* u = \w u+v -h u+ T^ u {t>0 -y^ 

that is uniformly bounded on each finite t-interval and satisfying the initial condition 
U(y, +0) = 1. It also satisfies the boundary condition U(+oo,t) = 1. In view of this, 
we consider a stationary solution S(y) = S^f*) that satisfies £(+00) = 1. For y > a, on 
the one hand, it is given by 



S(y, v *) = \f2irv*ye v * v I„{v*y) + 6K u (v*y) 



for some constant 9 G R. On the other hand we have two independent solutions e' 



and e a ~( y ~ a ) on (—00, a] (for v* > y/2/3/a), where a± = — v* ± Jv* — 2/3/a 2 , so that for 
some constants A + and v4_, 

%;0 = A + e a+{y - a) + A_e a - {y ~ a) . 

For the present purpose we have only to consider, as it turns out shortly, a solution with 
= 0, for which the continuity of S(y) and S'(y) at the joint a enforces 

_ a^S(a+) - S'(a+) -a + S(a+) + S'(a+) 
/i-L. — and J\— — . 



a_ — a + ct_ — 

We have the asymptotic formula 



y/2Trze- z I„(z) = 1 + fiz' 1 + 2~ x (l + /3)(3z~ 2 + 0(z~ 3 ) (z +00); 

we need to have the asymptotic form of the derivative, for which however we may simply 
differentiate term- wise ([2], p. 21). Hence, uniformly for y > a, as v* — > 00 

S(y) = l + ^- + 0{{yv*y 2 ) and S\y) = - J- - (1 ± ® P + 0{v^). 

yv* y 2 v* y A vi 

From this as well as a + = —f3/a 2 v* + 0(1/ v 3 ) we can readily infer that both of A± are 
positive and in particular S(y; v *) > 1 for all y, provided that is large enough. By a simple 
comparison argument we conclude that U(y,t) < S(y); in particular, U(y,t) = 1 + 0(1 /v) 
as desired. The proof of Proposition [13] is finished. 

Remark 4. The case x/t — > 00, < v < 1/2 can also be effectively analyzed by a formula 
obtained in [5], which for a = 1 may read 

* ( \a f'x-l ( (x-lf\ 

q ^ s)ds = JoT^^y — 2-^r 



, exp - ^ '- Us / -^e^^WV^ dy 

V2^? V 2s 1 1 



where 

I v (xy)K v (y) - I v (y)K v (xy) r°° -tx-i) y d V 



[K„(y)Y + * 2 'Mi-M A y' 

(In [5] this is proved also for 1/2 < v < 3/2 and similar equalities are obtained for the other 
values of v e R.) 
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